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Motivated by the presence of different orders in multilayered high-temperature super- 
conductors, we examine a model consisting of nonequivalent two Hubbard chains coupled 
by interchain hopping by using the density-matrix renormalization group (DMRG) and a 
mean-field theory. As an example, we consider a system with noninteracting chain without 
order and a Ifubbard chain with strong antiferromagnetic (AF) spin-density-wave correla- 
tion. We find that the magnitude of the interchain hopping controls the strength of induced 
AF correlation as well as that of original one. It is also found that the induced AF correla- 
tion decreases with increasing the magnitude of the original correlation. Implications to the 
multilayered system are discussed. 

§1. Introduction 

Recently the coexistence of antiferromagnetism and superconductivity lias been 
reported in multilayered high- Tc cuprate superconductors from nuclear magnetic 
resonance experiments.n''El' In the superconductors, there are different layers whose 
carrier concentrations are unequal. When the concentration is small (large), antifer- 
romagnetic (superconducting) order is dominating. At intermediate concentration 
region, the coexistence of both the orders is realized. The coexistence itself has the- 
oretically been examined by using models with strong correlation like a t-J model 
within an isolated single layer. ^' However, in addition to the possible coexistence 
in the single layer, we need to clarify a role of induced order from neighboring or- 
dered layers due to proximity effect, when we analyze the experimental results for 
the coexistence. 

In order to exploit the induced order, we consider a simple model that can 
describe the proximity effect: a two-leg ladder system with different interaction pa- 
rameters in each leg and coupled with one-particle electron hopping interaction. In 
this paper, we examine a two-leg ladder consisting of a noninteracting chain and 
of a Hubbard chain where antiferromagnetic (AF) spin correlation is dominating. 
It is noted that, though two-leg ladder systems with equivalent chains have been 
studied intensively by both analytical^ and numerical-- approaches, nonequivalent 
ladder systems have not been investigated as far as we know. For this purpose, we 
employ density-matrix renormalization group (DMRG) technique to calculate the 
distance dependence of spin-spin correlation function. We find that the magnitude 
of the interchain hopping controls the strength of induced AF spin correlation in the 
noninteracting chain as well as that of original correlation in the Hubbard chain. 
Furthermore, it is found that the induced AF correlation decreases with increas- 
ing the magnitude of the original correlation. This is reproduced by a mean-field 



E-mail: ydsumi@yukawa.kyoto-u.ac.jp 



typeset using VTpJ^.ds (Ver.0.9) 



2 



H. Yoshizumi, T. Tohyama and T. Morinari 



Ui, /ii 



leg-1 
leg-2 



Fig. 1. Schematic view of a two-leg Hubbard ladder. The hopping parameter t is common to both 
the leg 1 and leg 2, but the on-site Coulomb interaction Ui and the chemical potential fj,i depend 
on the l-th lee. 



treatment for the ladder system, indicating that these findings are independent of 
dimensionality. 

This paper is organized as follows. In §2, we introduce a model Hamiltonian 
of a two-leg ladder consisting of nonequivalent chains described by the Hubbard 
model and coupled with interchain hopping interaction. In §3, we present DMRG 
results and discuss how the interchain hopping affects on original and induced AF 
spin correlations on each chain. In §4, using an analytical calculation, we clarify 
a relation between the strengths of induced and original orders. §5 is devoted to 
summary of this work. 

§2. Model 

We investigate a ladder system consisting of two nonequivalent Hubbard chains 
connected by a hopping term, as schematically shown in Fig. [TJ The Hamiltonian of 
the two-leg Hubbard ladder model is given by 

n=Y,ni + n^ (2-1) 

1=1,2 

with 

i,cr i i 

and 

= -tp Yl (4,i,.c.,2,. + h.c.) , (2-3) 

1,(7 

where q ^ ^ i'^li a) annihilates (creates) an electron of spin a at the i-th rung on leg 

^ = Z^cr '>T'i,i,a = J2cr '^li a^i I ui ^ '^^ iutra-chaiu hopping parameter common to 
the legs, tp is the inter-chain hopping parameter, Ui is the on-site Coulomb repulsion 
on leg /, and is the chemical potential that controls the charge density on each 
leg. Hereafter we take t = \. 

When tp = 0, the ladder represents decoupled chains. The parameters in the 
chains determine dominant correlations among spin, charge, and superconducting 
orders at half filling.''' In this study, we consider a simple case where the I = 1 leg 
is noninteracting {Ui = 0) and the I = 2 leg has positive number of U2 giving rise 
to the longest power-low correlation for spin channel. In the following, we call the 
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I = 1 leg the metal chain and the / = 2 leg the spin-density-wave (SDW) chain. We 
connect the two chains by introducing tp, whose value is assumed to be small, i.e., 
tp < t, keeping the multilayer cuprates in mind. In order to escape the change of 
correlation functions due to charge transfer from one chain to another, the chemical 
potential m is chosen to make each chain contain a half of total electrons, i.e., the 
averaged electron number in the I leg is given by ^i{ni^i) = N/2, where is the 
number of total lattice sites and (• • • ) stands for the ground-state expectation value. 
In this study, we consider the case that the total number of electrons is equal to A^, 
i.e., half filling, and take ni = Ui/2. 

Employing the finite-size algorithm of density matrix renormalization group 
(DMRG) method^ for a ladder with = 64 x 2 sites under open boundary condi- 
tions, we calculate the spin-spin correlation function for each leg, which is defined 
as 

Q(|r,-r,|) = (5f(r,)Sf(r,)>, (2-4) 

where Si{ri) is the z component of the spin operator at the i-th position rj on the 
leg I. The middle of the site i and the site j is chosen to be at the center of the leg 
to minimize size effects coming from edges of the ladder. 

§3. DMRG results 

3.1. Spin-spin correlation in the SDW chain 

We first examine spin-spin correlation in the SDW chain C2{r) (the I = 2 leg). 
Figure [2] shows the correlation function for several values of the interchain hopping 
tp in the case of U2 = 4. We plot only positive values of the correlation where the 
distance r is of even number. We note that odd distances give negative values of the 
correlation. 

Before discussing the details of the results, we should mention convergency of 
the data in terms of the truncation number m in DMRG. Since we treat a ladder 
system with a noninteracting chain and small tp, we are required to use a large 
number of m. Inset in Fig. [2] shows m dependence of C2(r) for the case of tp = 0. 
Since the two chains are decoupled, C2(r) must agree with that of a 64-site single 
chain with on-site Coulomb interaction U = 4, which shows a linear dependence in 
the log- log plot for r < 40 (solid squares). With increasing m from 400 to 800, C2(r) 
approaches that of the single chain. The data at short distances (r < 10) seems to 
converge in the logarithmic scale for m = 800, while the data at long distances show 
slow convergencey. Such slow convergency comes from the small correlation between 
the two chains. The convergency in terms of m becomes much better with increasing 
tp. However, m = 800 are not enough to get good convergency at large distances. 
Unfortunately m = 800 is a maximum number we can access using our available 
computers. We note that CPU time for the calculation of the correlation function 
for a given parameter is about 52 hours by using 64 cores on FUJITSU HX600. We 
also note that truncation error for m = 800 is about 10~^. 

Keeping in mind the convergency problem mentioned above, let us discuss how 
C2(r) is affected by the introduction of tp. For r < 10, C2(r) slightly decreases with 
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Fig. 2. Spin-spin correlation function on the SDW chain C2{r) as a function of the distance r 
between two sites for a 64x2 ladder with U2 ~ 4. The data for several values of the interchain 
hopping tp are shown. The truncation number m = 800. Inset: m dependence of C2{r) for 
tp — Ohy using the ladder system. The filled squares represent the case of a 64-site single chain. 



increasing tp. In order to clarify the effect of the on-site Coulomb interaction, we 
show in Fig. tp dependence of the spin-spin correlation at r = 4 on the SDW 
chain for U2 = 2, 4, and 6. For the check of convergency with respect to m, we 
show 1/m dependence of C2{r = 4) for C/2 = 4 in the inset. The data for tp = 
seems to be safely extrapolated to the data of a single chain plotted at 1/m = 0. 
Similarly, the data for tp = 0.4 exhibits a smooth change toward 1/m = 0. These 
indicate good convergency for m = 800. In the main panel of Fig. [31 C2{r = 4) is 
suppressed with increasing tp. This is easily understood as a result of mixing with the 
connected chain. Such suppression has been discussed in the context of the effect of 
apical oxygen on pairing amplitude in the Cu02 plane,!^' though the dimensionality 
is higher than the present case. The suppression of C2{r = 4) depends on U2- With 
increasing U2, the suppression is reduced and very small at U2 = 6. This means 
that, at large U2, C2 hardly changes even if tp is introduced. This is also easily 
understood by the fact that the double occupation on the SDW chain is not allowed 
during the process of interchain hopping. This effect is also seen in the metal chain 
as a suppression of induced AF correlation, as will be discussed below. 

At large distances more than r ~ 10, C2{r) increases with increasing tp as 
shown in Fig. O Since the convergency of the data at this region is questionable 
as mentioned above, we cannot give a definite answer to a question whether C2{r) 
increase with tp or not. We need larger m in DMRG or analytical studies in order to 
clarify this question. This remains as a future work. We note that a linear behavior 
up to r ~ 30 in the log-log plot for the case of tp = 0.4 indicates no spin gap in the 
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Fig. 3. tp dependence of spin-spin correlation on the SDW chain at a selected distance r — 4 for 
a 64x2 ladder and different values of U2- The truncation number m = 800. Inset: C2{r — 4) 
for ?72 = 4 as a function of the inverse of m. The data for tp = and 0.4 are shown. The filled 
square at 1/m = is the data obtained from a 64-site single chain. 



SDW chain. However, we need further study to draw a final conclusion. 

3.2. Spin-spin correlation in the metal chain 

We show in Fig. [5] induced AF spin correlation function ACi in the metal chain 
(the I = 1 leg) , which is defined as the difference of the correlations for finite tp and 
tp = 0. The data at even distances are plotted. At U2 = 4, the magnitude and length 
scale of ACi increase with increasing tp as shown in Fig. H^a) . As for convergency of 
the data in terms of m, we again plot 1/m dependence of ACi at several values of 
r for tp = 0.3 in the inset of Fig. IHa). The data gradually converges from m = 500 
(1/m = 0.002) to m = 800 (1/m = 0.00125). In Fig.Hl^b), we plot U2 dependence of 
the induced AF correlation for a fixed value of tp = 0.3. With increasing U2, ACi 
increases and saturates near U2 = 3. Then, ACi decreases with further increasing 
U2. 

In a single Hubbard chain, the magnitude of the AF spin correlation increases 
with increasing U. The increse of U also reduces charge fluctuation. This relation 
between AF correlation and charge fluctuation tempts us to relate the nonmonotonic 
U2 dependence of the induced AF correlation ACi to charge fluctuation in the metal 
chain. In order to confirm this, we have calculated the charge fluctuation on a central 
site of the metal chain, defined as (An)'^ = {nf^^2i=-\) ~ {^1=32, 1=1)'^ (not shown). 
As expected, (Z\n)^ decreases from = 2 to 3, and then increases with further 
increasing U2, although the change of its magnitude is only 0.5%. This implies 
that the induced AF spin correlation is related to charge fluctuation affected by the 
connected SDW chain. 

In order to make more clear the U2 dependence of ACi, we examine a relation 
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Fig. 4. Induced AF spin correlation function on the metal chain ACi (r) as a function of the distance 
r between two sites for a 64x2 ladder. The truncation number m — 800. (a) U2 ~ 4 and the 
data for several values of tp are shown. Inset: ACi{r) at several distances for tp — 0.3 as a 
function of 1/m. (b) tp = 0.3 and the data for several values of U2 are shown. 



between the original AF spin correlation on the SDW chain at tp = and the induced 
correlation at tp = 0.3 on the metal chain. For this purpose, we define the following 
integrated correlation functions: 



4'^= E C2{\ri-rj\), 



i~j=evcn 



i—j=even 



(3-1) 
(3-2) 



where i—j runs over even numbers, i.e., i and j are on the same sublattice. Therefore, 



r(2) 



is roughly proportional to the magnitude of spin structure factor at momentum 

r(l) 



Q = vr on the SDW chain. Similarly, Alfj roughly corresponds to the induced 

part of the Q = it structure factor on the metal chain. We show 'AI^^ versus I^^ in 
Fig. O Each data point corresponds to a given value of C/2. We find that the induced 
AF correlation shows non-monotonic dependence on the original AF correlation on 
the SDW chain. In particular, the induced correlation has a maximum for a certain 
value of the original correlation and weakens with further increasing the strength of 
the original correlation. This behavior seems to be counterintuitive, since induced 
correlation is expected to be strong if original correlation is strong. In the next 
section, we consider the origin of this behavior using a mean-field treatment of the 
ladder system. 



§4. Mean-field analysis 



We perform a mean-field SDW approximation for the SDW chain (the I = 2 leg). 
By assuming an ordered momentum Q = n under periodic boundary condition, 
the corresponding SDW order parameter is given by A2 = ^(nj2t — ni2j,)(— 
which becomes independent of site i. Using A2, we obtain the Hamiltonian of the 
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Fig. 5. Relation between integrated AF spin correlation I^' 

integrated induced AF spin correlation Alj-V on metal chain at tp = 0.3 for a 64x2 ladder. The 
data are taken from DMRG results with m = 800. 



nonequivalent two-leg Hubbard ladder: 

'Hmf = + ^2 + 'H± 



with 



^2 = X;4,2 



/ 6 \ 

0^ 



\ Ck+Q j 

I ik -C/2Z\2 \ 

Cfc C/2Z\2 

ik+Q 

V U2A2 ) 



0'k,2, 



and 



fc.lT 



where ^k = — 2t cos k—^ with the chemical potential ^ and ^ = (c^ ; -|-, ; j.' "^fc+Q i f ""fc+Q ' i 
with ^ ^ being the Fourier component of cj^ ^. Integrating out the freedom of the 
SDW chain, the inverse of an effective Green function on the I = 1 leg reads 



(4-1) 



(4-2) 

(4-3) 
(4-4) 

(4-5) 



/ Bn{k,A2) -Dn{k,A2) \ 

Bn{k,A2) Dn{k,A2) 

-Dn{k,A2) Cn{k,A2) 

V Dn{k,A2) C„(A;,Z\2) / 



, (4-6) 
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where Bn{k, A2) = Ck+An(k, A2) {itOn - ik+q), Cnik, A2) = Ck+Q+An{k, A2) (iw„ - £,k) 
Dnik, A2) = U2A2An{k, A2), and Anik, A2) = tl/[{iu;n-^k)iii^n-^k+Q)-iU2A2f] , 
with LOn being the Matsubara frequency 

By relating the off-diagonal elements in (j4-6p to the induced order, an effective 
SDW order emerging on the / = 1 leg, Af, is given by Dn{k, A2) /U2- Considering 
electrons near the Fermi level at half- filing, i.e.. A; = 7r/2 and /Li = 0, and assuming 
tp <C Z\2, we obtain 

where the induced order decreases inversely proportional to the magnitude of the 
original order. This behavior is qualitatively consistent with the DMRG results of the 
induced AF spin correlation as shown in Fig.[5j Of course, the quantities calculated in 
the DMRG and mean- field calculations are different, but underlying physical origins 
should be the same because both the correlation and the order are related to the 
interplay of on-site Coulomb interaction and the spin degrees of freedom. From the 
mean- field result, we can obtain an intuitive picture of the decrease of the induced 
order at large original order. For large original order, an electron in the metal chain 
feels strong potential due to the order and thus the electron is hard to hop to the 
SDW chain. As a result, the transfer of the order to the metal chain is less effective 
during the second order process using tp as seen in (|4-7p . It is interesting to note 
that this mechanism is independent of dimensionality. 



§5. Summary 



We have examined the spin-spin correlation function of a nonequivalent two-leg 
Hubbard ladder in order to exploit the induced order due to the proximity effect. 
One leg is a noninteracting chain and the other is a Hubbard chain with positive on- 
site Coulomb interaction where AF spin correlation is dominating. The two chains 
are coupled by interchain hopping. By employing the DMRG technique, we find 
that the AF correlation on the SDW chain is slightly suppressed at short distances 
with increasing interchain hopping. This suppression is due to the transfer of AF 
spin correlation to the metal chain. Then, the induced AF spin correlation in the 
metal chain is enhanced. Fixing tp, we find that the induced correlation decreases 
with increasing the magnitude of the original correlation in the SDW chain. This is 
supported by a mean-field SDW calculation. 

The model taken in this study is too simple to describe multilayered high- Tc 
cuprates in which the coexistence of antiferromagnetic and superconducting orders 
has been reported. The issue of the coexistence is beyond the present study, but a 
possible induced order due to proximity effect can be discussed based on the present 
results of DMRG and mean-field calculations. The results indicate that, if the order 
is large in one Cu02 plane, the induced order in neighboring planes should be small. 
This means that, if the inner Cu02 plane in the multilayered cuprates is close to half 
filling where the magnetic moment is large, induced part of the moment in the outer 
plane should be small. The present results might contribute to further understanding 
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the coexistence of different orders in the multilayered cuprates. 
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